Composite fermions in the Fractional Quantum Hall Effect: Transport at finite 

wave vector. 
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We consider the conductivity tensor for composite fermions 
in a close to half-filled Landau band in the temperature 
regime where the scattering off the potential and the trapped 
gauge field of random impurities dominates. The Boltzmann 
equation approach is employed to calculate the quasiclassical 
transport properties at finite effective magnetic field, wavevec- 
tor and frequency. We present an exact solution of the kinetic 
equation for all parameter regimes. Our results allow a con- 
sistent description of recently observed surface acoustic wave 
resonances and other findings. 

PACS numbers: 71.10.Pm, 73.50.Bk, 73.20.Dx 



The properties of two-dimensional (2D) electron sys- 
tems in high magnetic fields appear to be well described 
by the model of composite fermions (CF's) Q]^]. It pos- 
tulates the existence of fermionic quasiparticles consist- 
ing of electrons to which an even number of flux quanta 
has been attached. The corresponding mapping can be 
expressed by introduction of a Chern-Simons (CS) gauge 
field The average magnetic field B q q experienced 

by the CF's is the difference of the external magnetic 
field B and the mean field of the flux tubes, = (fxfion, 
where 4>q — hc/e is the flux quantum, <f> is the (even) 
number of flux quanta per electron and n is the elec- 
tron density. The CF model appears to capture the ef- 
fect of the Coulomb interaction in maintaining optimum 
electron separation particularly near even denominator 
fillings of the Landau level, v = l/(2q), when choosing 
<j> = 2q (q = 1, 2, . . .). Then the effective magnetic field 
B„q is largely reduced near and exactly zero at these fill- 
ing factors. In this regime the CF model implies the exis- 
tence of a Fermi surface with wave number kp — y 'Aim in 
a weak magnetic field ||. At low temperatures the CF's 
are subject to elastic impurity scattering off the random 
potential created by a remote layer of donors and, more 
importantly, off the frozen random gauge field configu- 
ration induced by the impurities. In addition, the CF's 
interact via the scalar (Coulomb) potential and the CS 
gauge field. 

The existence of new particles with a definite effective 
mass moving under the influence of an effective mag- 
netic field B c q has been demonstrated in a number of 
experiments. For example, the interpretation of magne- 
tooscillations in the longitudinal resistivity in terms of 
the Shubnikov-de Haas (SdH) effect for CF's has been 
very successful H,p| . The observation of resonance effects 



for the motion of charge carriers in the field of surface 
acoustic waves (SAW) is another example ||. However, 
a closer examination suggested [jfl that the experimental 
observation of these resonances is inconsistent with the 
conventional interpretation in the range of wave frequen- 
cies used. A quantity of central interest for the theo- 
retical description of this phenomenon is the conductiv- 
ity tensor at finite wave vector a^lu, q), the calculation 
of which we will reexamine in the present paper. The 
momentum-dependent conductivity a ^(us , q) , which de- 
termines the gauge field polarization operator, is also of 
importance for the damping of SdH oscillations by static 
and dynamic fluctuations of the gauge field |||| • 

In this letter we present a quasiclassical calculation 
of o {to, q) in a magnetic field, taking into account the 
impurity scattering of fermions off a static random mag- 
netic field and random potential. Surprisingly, despite 
the general character and the importance of the prob- 
lem, we did not find a correct calculation of er M „(w,g) 
in the literature even for the usual case of a short-range 
correlated random potential, except in the collisionlcss 
limit fl. The standard reference jio) appears to have 
missed the fact that, since at finite wave vector q the 
current and density fluctuations are coupled, the relax- 
ation time approximation in its simplest form cannot be 
used because it violates the particle number conserva- 
tion. In particular, the results of Ref. (used later for 
interpretation of the SAW experiments in Refs. |§,^|,0) 
violate the exact property lim^^o q^cr ^(u: , q) = 0, which 
is a consequence of the current conservation. The par- 
ticle number conservation has to be taken into account 
properly in modelling the collision integral of the Boltz- 
mann equation, as shown below. Another new aspect of 
the transport problem we are considering is the extreme 
angular dependence of the scattering off a random mag- 
netic field, which is so strongly peaked in the forward 
direction, that the relaxation time approximation is not 
adequate. 

In quasiclassical transport theory the components of 
the electrical current may be expressed in terms of 
the distribution function f(k;r,t) of fermions of mo- 
mentum k at position r and time t as j^{r,t) 
• f(dk)(k^/m*)f{k;r,t), where e and m* are the charge 
and the effective mass of the carriers. We are interested in 
the linear response fi(k;r,t) to an applied electric field 
E{r,t) = i£exp(— iiot + iqr) in the presence of a weak 
magnetic field Bo perpendicular to the plane, such that 
lo c t s <C f (here lo c = eBo/m*c is cyclotron frequency and 
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t s 1 is the width of the Landau levels). The function /i 
is obtained from the kinetic equation 

i(w - vq)h + ~{v x B ) • Vfc/i + e^V fc /o = C{/i}, 

c 

where /o is the equilibrium distribution function, v — 
k/m* (we assume isotropic Fermi surface) and C{f%} is 
the collision integral. The external force oc Ev(df/de) 
acts only on states at the Fermi circle (e = ep), and 
therefore fi is only a function of the angle <j> of k with 
the x axis, say. In terms of the angular momentum eigen- 
functions Xm(^) = e lm ^, the collision integral is given as 

C ifi}= E A »> / 7rexp{*m(0-0')}/i(A (1) 

Particle number conservation requires Ao = 0, also A m = 
A— m = A,^, since C is a hermitean and non- negative 
operator. The usual relaxation time approximation used 
in Ref. pTcj ), C ~ /i/t, is equivalent to setting A m = 1/r 
for all m, obviously violating the requirement Ao = 0. 
This can be corrected by approximating 

C~(l/r)[/i- //i#/2vt], (2) 

implying Ao = and A m = = 1/r for all |m| > 1 (in the 
following called model I). Model I describes isotropic im- 
purity scattering and is applicable for short-range cor- 
related random potential. For scattering off a random 
magnetic field, when the transition rate W(4> — </>') for 
scattering of a particle from angle </)' to angle (j> on the 
Fermi circle is given (in Born approximation) by Jll]] 
W((t) - <f>') = r~ 1 cot 2 [(0 - 0')/2], model I is not ade- 
quate. In this case, W is seen to diverge in the forward 
scattering limit <f> = (p 1 . The eigenvalues of the collision 
operator 

C{/i} = / (d4>'/2n)W(4> - 4>')[h (<£) - h(<t>')] (3) 

are then given for m ^ by A m = T -1 (2|m| — 1), i.e. they 
increase linearly with \m\ (to be called model II). The 
larger \m\, the deeper the forward scattering divergence 
is probed. Finally, if both random potential and ran- 
dom magnetic field scattering mechanisms are present, 
we have A m = t -1 (1 — 2p + 2p\m\), with < p < 1, 
where p = and p — 1 correspond to pure models I and 
II respectively. In the sequel, we solve the Boltzmann 
equation exactly for arbitrary values of the frequency ui, 
momentum q and magnetic field Bo, and for arbitrary 
value of p. 

Solution of the kinetic equation. We define a di- 
mensionless distribution function g(<p) by fi(k;q,u>) = 
(—dfo/de)elEg((j)), where I = vt is the transport 
mean free path. The function g((f>) obeys the integro- 
diffcrential equation 

ir(uj — vq cos <j>)g — uj c t— — cos(0£ — (f>) = rC{g}, (4) 



where tan</>£ = E y /E x and we have chosen the direc- 
tion of the momentum q || x. Expanding g in terms of 
eigenfunctions Xm(<l>), g{4>) = Em=-oo5™Xm(^)> yields 
the recursion relations 

a m g m - b{g m+ i + g m -i) = S m (5) 

with a m = —iujT + imuj c r + A m r, b = —ivqr/2, and 
S m = (1/2) exp(— im<j>E)S\ m \ i- We note that the "an- 
gular momentum" components g m only mix through the 
term proportional to q. For positive m > 2, the system 
(0) can be solved to yield 

g 2 = R +9l ■ R+ = (b/a 2 )P 3 /P2 , (6) 
where Pj~ are power series in b 2 (thus, in q 2 ): 

oo n-\-k— 1 

P fc = ^(l/n!)(-6 2 /a + )7 II a < ( 7 ) 

n=o i=k 

Here a+ = \\m m ^ co {a m / m) , and we assume a m to in- 
crease at most linearly with m as m — > 00; the con- 
dition which is fulfilled for the models considered. For 
negative m we find analogously <?_2 = R—g—i, where 
i2_(w c ) = R+{— lo c ) and a-(uj c ) = a + (— uj c ). In terms 
of the components g\, g-\ of the distribution function 
the components of the conductivity tensor are given by 
o-xn = 00(51 +5-1) and &vn = icr (cji - g~i), where 
(To = e 2 nr /m* and the index /i corresponds to the direc- 
tion of the electric field determined by the angle 4>e- Us- 
ing eq.(Q) and its counterpart for negative to, we solve the 
system of the remaining three equations from (^) (with 
to = 0, ±1) and find 

<7xx = <7 kjt(7 + + j-)/D ; 

a xy = <7 u/t(7_ - l+)/ D = -&yx ; (8) 

(T yv = a [4b 2 + ico T (j + + j-)]/D , 

where D/2 = 6 2 (7+ +7-) + iu>rj + j^ and y± = 1 — iuTzL 
iu c T — bR±. For the chosen direction of q \\ x we find 
as expected a xf i(q,i^) oc uj for lo — > 0, q ^ 0. The exact 
solution (^|) of the Boltzmann equation for any impurity 
collision integral can be easily evaluated for any finite 
value of q, since the power series representation of the 
functions P/., eq.(0), is rapidly converging. For the spe- 
cial cases of both models I and II the eigenvalues of the 
collision integral A m and consequently the coefficients a m 
are linear function of to, so that the power series repre- 
sentation of Pk is reduced to that of the Bessel function 
J v (z), with complex valued index v and argument z: 

R± = J2+p±(Q±)/Jt+p ± (Q±); Q± = -ivqr/a±, (9) 

with a± = 2p ± iuj c T, j3± = (1 — 2p — iujr)/a±. 

Substituting eq.(^|) in (g) and using the recursion re- 
lation for the Bessel functions, z[J u+ i(z) + J u -i{z)\ = 
2vJ v (z), we finally get the conductivity tensor for arbi- 
trary values of the three dimensionless parameters Q = 
qR c = qv/uj Cl S = uj c t, and T = uj/uj c : 
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(ia /QS)(J-J 1+ - J+Jj-) 
J-J+ - (Q/2T)(J_J 1+ - J+ J x _) ' 
(a a /QS)(J^J 1+ + J+J^ 



>yx 



'yy 



'xy 



J-J+ - (Q/2T)(J_ Ji+ - J+ Ji. 
(2i f 7 /5T)J 1 _J 1+ 



J- J+ - (Q/2T)(J_ Ji+ - J+Ji_) ' 
J± = •^±[T+(i-2 P )i/s]/(i±2p l /s)( ( 5/(l ± 2pi/S)): 
Ji± = Ji±[T+(i-2 P )i/s]/(i±2 P i/s) (Q/ (1 ± 2pi/ S)) (10) 

In the limiting case p = (model I) the formulas ([it)]) sim- 
plify, since the arguments of all the Bessel functions are 

equal to Q: J± = J±(T+i/s){Q), J i± = J i±(T+i/s)(Q)- 
In this case eqs.(|10|) can be rewritten in a slightly modi- 
fied form (similar to that obtained in Ref. Q in the limit 
t = oo) by using the Bessel functions identity 



J- J 



i+ 



J+Ji- 



2(T + i/S) 



J-J 4 



sin7r(T +i/S) 
tt(T + i/S) 



For the particular case of model I, the kinetic equation 
([|), (Q) can be also solved in a different way. With the 
notation f (d4>/2Tr)g(4>) — go, eq.(|J) reduces to a first 
order differential equation, which has the solution 



/ d^'exp^OM')} 



.go - cos(0 B 



(11) 



ivq 



(sin0-sin<//) (12) 



Calculating now J d<fig(<fi), we get a selfconsistency equa- 
tion for go. Substituting the found value of go back into 
cq.(|ll]), and calculating the conductivity components, we 
find 



LOT I AqO 

a xx = -2i<r - ~ I 1 - iur 



Pyx — ~0~xy — —2<7q 



(qvr) 2 \~ A 00 - lo c t 

LO A s0 



vq A o - lo c t 
2ao ( AosA s o . 
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lo c t \Aoa - lo c t 

2tt 



(13) 



Aij = 7T / rf ^(0) / #Vi(0')ex P {X(0,0')} ; 

27r JO J-oo 

4>o(</)) = l ; ips(<P) = sin0 

In this form the results for model I are more suitable for 
consideration of the zero magnetic field limit (see below) . 

Limiting cases. Despite their apparent simplicity, the 
results (pi) possess a reach variety of regimes with vari- 
ous types of behavior depending on the relations between 
S, Q, and T; we consider only some of them below. In 
the limit q = we have R± — 0, so that the general 
results (||) reduce to the usual Drude formulas. In the 
collisionless limit r = oo, formulas (|o|) can be proven to 
reduce to the results of Simon and Halperin ||. 

In zero magnetic field, lo c = 0, we have o~ xy = and 



2ia. 



P {i-i P )Wi + ^/(i-i P y 



yy 



JI 



1 - (1 - ip)y/l + K 2 /(l-l P ) 2 
2aoWl + K 2 /(l - ipf - 1)(1 - i P )/k 2 



2op 

K 



(_ip-l)/2 



(«) 



f(l-»p)/2(K) 
2<7q I(l-ip)/2( K ) 
« ^(-ip-l)/2(«) 



K 
IP 



(14) 



where the superscripts / and II refer to the models I 
and II respectively, K = qvr, p — lot, and I v {z) is the 
modified Bessel function jl2|]. In particular, in the large 
wavevector limit, k» l,p, we get for both the models 
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-ie vLo/q ; 



e 2 w/q , 



(15) 



where v = m* /2n is the density of states at Fermi sur- 
face. These results are identical to those in the collision- 
less limit. We conclude that although the eigenvalues A m 
of the collision operator C of model II grow with m with- 
out bound, for qv ^> t~ 1 ,lo c the collision integral C is 
negligible in comparison with vq. 

Finally, let us consider the case of relatively strong 
magnetic fields, S = lo c t 3> 1, corresponding to the re- 
gion where most of the resonance and magnetooscillation 
effects for the composite fermions are observed JllJ]. For 
these phenomena the behavior of conductivity at low fre- 
quencies to <C Lo c ,qv and finite wavevectors qR c > 1 is 
important. In this regime, a xx and a xy are small (since 
they vanish at to = 0), and a yy plays the most important 
role. We find 



a yy ~ 2a Jt(Q)/(l - J'oiQ) - iwr) 



mi 



<?oJi(Q)/Q 2 



J 2 (Q) + J Q 2 (Q) 
2a /QS , 



MQ)MQ)/Q ' 



Q/S<1 
Q/S^i 



Note that for the model II at a large wavevector, Q ^> 
S, the zero magnetic field result, eq.(15), is restored, 
whereas it is not the case for the model I. 

SA W resonances. As an important example of appli- 
cation of our results, we will now reconsider the interpre- 
tation of experimental data on the SAW resonances [pUTf- 
As discussed in Ref. {jj , using the formulas from Ref. |ly] 
for the CF conductivity does not allow to explain the 
observation of resonances at the experimental frequency 
lo = 2tt x 10.7 GHz if the value of the effective mass 
extracted from the SdH measurements, m* ~ 0.8m e , is 
assumed (m e being the free electron mass). We have, 
however, recalculated the SAW velocity shift nM, 



Av/v = (a 2 /2)Re[l + ia^/a-n 
Vxx - Pyy/Ply Pxy = 2h/e 2 ; 

Pyy = a xx / {o~xx°~yy + a X y) 



p 1 ; a 2 j2 = 3.2 • 10" 4 ; 

(16) 



making use of the correct formulas for a ^(ui , q) derived 
above and setting m* = 0.8m e . The constant a m in 



3 



eq. (|T^) and the transport time r were used as parameters 
to optimize the fit to the experimental data from Fig.l 
of Ref . . We get very good fits to the experimentally 
observed resonances by using any of the models I or II, or 
a mixed model with combined mechanism of scattering 
(Fig.l), at a rn = 1 • 10 6 cm/s ~ Kr 6 rr\ which is close 
to expected values of this parameter . Note that the 
value of the transport time r found from optimizing the 
fit depends appreciably upon the model assumed, varying 
from t = 20 ps to r = 120 ps in the limiting cases of the 
models I and II respectively. 



> 
< 




FIG. 1. Velocity shift of the surface acoustic waves, eq.(|l6|), 
as a function of the effective magnetic field experienced by 
composite fermions, with to the following values of the param- 
eters: to = 2tt-10.7 GHz, q = 2tt/(2.7- 1CT 5 cm), m* = 0.8m e , 
n = 1.6 ■ 10 11 cm -2 , and a m = 1 ■ 10 6 cm/s. The three curves 
shown correspond to: a) model I (p = 0), r = 20 ps; b) mixed 
model (p = 0.3), r = 50 ps; c) model II (p = 1), r = 120 ps. 
For clarity, the curves a) and b) are off-set by 0.0001 and 
0.00005 respectively. The principal and secondary resonances 
(minima) are seen at Bgfj ~ 6 kG and -B e ff — 2.5 kG, as 
observed in the experiment JtJ. 

In conclusion, we have presented a solution of a very 
general problem: calculation of the scmiclassical conduc- 
tivity tensor of fermions in a random potential and/or 
random magnetic field in 2D at finite values of frequency 
lu and momentum q and in uniform magnetic field Bq. 
We have pointed out that explicit account of the particle 
number conservation in modelling the collision integral 
is crucially important at finite q; missing this fact led to 
incorrect results in earlier publications. Expanding the 
distribution function in angular harmonics in momen- 
tum space, we mapped the Boltzmann equation onto a 
system of recursion relations and found its solution. In 
the particular case of a short-range correlated random 
potential and/or random magnetic field the conductivity 
tensor takes a rather simple form, eq.(10). Application 
of our results resolves an apparent inconsistency with the 
interpretation of the experimental data on resonances of 
composite fermions with surface acoustic waves. Using 



our formulas with experimental values of frequency and 
the CF effective mass m* — 0.8m e , as extracted from 
the SdH measurements, we obtained the curves (Fig.l) 
which describe well the experimentally observed reso- 
nances (Fig.l of Ref. 0). Our results are also to be 
used in all the cases when the propagator of the gauge 
field (coupled to CF's) at finite uj and q is considered. 
This is important, in particular, when one discusses the 
damping of SdH oscillations at finite temperatures due 
to dynamic fluctuations of the gauge field ||U ■ 
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